Signature of Quantum Criticality in the Density Profiles of Cold Atom Systems 
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In recent years, there is considerable experimental effort using cold atoms to study strongly 
correlated many-body systems [T}{6] . One class of phenomena of particularly interests is quantum 
critical (QC) phenomena. While prevalent in many materials, these phenomena are notoriously 
difficult theoretical problems due to the vanishing of energy scales in QC region. So far, there are 
no systematic ways to deduce QC behavior of bulk systems from the data of trapped atomic gases. 
Here, we present a simple algorithm to use the experimental density profile to determine the T=0 
phase boundary of bulk systems, as well as the scaling functions in QC regime. We also present 
another scheme for removing finite size effects of the trap. We demonstrate the validity of our 
schemes using exactly soluble models. 



In a recent paper, we gave a general scheme to derive 
thermodynamic functions of bulk systems at finite tem- 
perature from the density profile of trapped gases within 
local density approximation (LDA)[7 . This scheme has 
recently been implemented by Salomon's group at ENS to 
map out the global phase diagram and equation of state 
of unitary Fermi gasj^^^. The study of quantum phase 
transitions, however, involves other challenges. First of 
all, these are phase transitions at T = caused by 
changes of the parameters of the system. In addition, 
there are systems with phase transitions only at T = 
but not at finite temperatures. How to access the T = 
phase diagram from the data of trapped gases at non-zero 
temperature is not obvious. Secondly, cold atom experi- 
ments typically have 10^ to 10^ particles, which will give 
rise finite size effects. How to eliminate these effects sys- 
tematically is also unclear. The purpose here is to offer 
solutions to these problems. 

(A) Algorithm for uncovering quantum criticality from 
non-uniform density profile: Let us consider the equation 
of state n = n(T^ ji)^ where n is the density and /i is the 
chemical potential. Let jic be the chemical potential at 
T = where the phase transition takes place. Near ji = 
/ic, general arguments show that when the dimensionality 
d is below the upper critical dimension dc^ n = n{T^fi) 
becomes [10] 

n(M,T) - n,(T,M) = Ti+'-^G , (1) 

where u is the correlation length exponent, is the reg- 
ular part of the density, and Q{x) is a universal function 
describing the singular part of the density near criticality. 
(Above dc^ d/ z -\- 1 — 1/uz is replaced by the mean field 
value 1/zu). Correspondingly, the compressibility = 

dn/dj^ obeys T-^-^ + ii [/^(/i, T) - /^^(/i, T)] = (^) 

for d < dc. For d > dc, — 1 — f + ^ is replaced by 0. 

Eq.Q is particularly useful when rir of one of the 
phases is known. For boson and femion Mott insula- 
tors, Ur are integers (hence /^^ = 0). For systems such as 



the fully spin polarized fermi gas, can be calculated 
analytically. In such cases, if we plot the "scaled den- 
sity" A(/i,T) = T-i-^/^+i/^^(n(/i,T)-n^), or the scaled 
compressibility C(/i,T) = + A [A>:(/i, T) - nr{iJi,T)] 

versus /i for different temperatures, then all curves will 
intersect at the same point ji = jic Once jic is deter- 
mined, one can then plot A{ii^T) (or C{ii^T)) versus 
u = {ji — \ic)lT^I^^ . The scaled density curves for all 
temperatures will collapse into a single curve, which is 
the scaling function Q(u) (or Q'iu)). 

This scheme can easily be implemented in cold atom 
experiments within LDA, which takes the measured den- 
sity rf-^ (x) at point x at temperature T as its bulk value 
n{/jj^T) with /i replaced by /i(x) = /i — F(x), where 
F(x) is the confining trap; i.e. n^^(x) = n(/i(x),T). 
With T and /i determined from the measured density 
profile[7| lUj, n = n^ji^T) of the bulk system is read- 
ily obtained from the non-uniform density profile n^^(x) 
of the trapped gas. To locate the critical point /ic, one 
simply plots A^^(r) = T-i-^/^+i/^^(n^^(r) - n/) ver- 
sus /i(r) for different samples at different temperatures 
T and looks for their common intersection. Such inter- 
section will only occurs if z and v assume their correct 
values. Recalling that z is typically either 1 or 2, and v 
can be estimated from our knowledge of various univer- 
sality classes, the number of trials needed in practice to 
produce the common intersection is very small. 

Even if rir is unknown, its regularity allows one to ex- 
pand it around (/i = jic^T = 0), nriT^ji) = n^(0,/ic) + 
Z)(n,m)^(o,o) ^n,mT'^{l^- l-ic)'^' One cau then implement 
the previous scheme by retaining the first few terms in 
the series, treating their coefficients an,m as parameters 
to be adjusted to obtain a crossing of the scaled densities 
or compressibilities at various temperatures. 

(B) Quantum criticality of ID hard core bosons: To 
demonstrate our scheme, we use ID hard core bosons and 
ID free Fermi gas as examples. Since both systems have 
exact solutions in confining traps and in free space, we 
can then characterize the accuracy of LDA by comparing 
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FIG. 1: Phase diagram of ID hard core Bose gas. / and // 
are the quantum critical points of vacuum (V) to Tomonaga- 
Luttinger Liquid (TLL), and Mott (M) to TLL transitions. 
The quantum critical regions associated with quantum critical 
point / and // are denoted as QCI and QCII, given by |/i + 
2J| < T and |/i — 2 J| < T respectively. The shaded region 
in the vacuum and the Mott phase represent the low fugacity 
regime of particle and holes. 



the density calculated from it with that calculated from 
the solution in a harmonic trap, (denoted as n^^(x)), 
which now plays the role of experimental data. Our 
hamiltonian is H = — X](r r') (^r^r' + ^'^') + "^R^R^ 
where R labels the lattice sites , (R, R') are nearest 
neighbors, 6^ creates a boson at R, V^i = ^Mcj^R^ 
is the harmonic potential, and nR = b^J:)-^. The hard- 
core constraint is implemented by restricting nR = 
or 1. It is well known that the wavefunction for ID 
hardcore bosons, ^ b{Ri-,R2-, '"•,Rn)i is that of a free 
Fermi gas R2-, - Rn) multiplied by a factor A = 

\Ir>r' ^{R-R'), where e{x) = 1 or (-1) if x > (x < 0); 
i.e. = A^p. The density profile of hard core bosons 
is then given by that of the free Fermi gas. 
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\u^{x)\^f{En), 



(2) 



where Un is the eigenfunction of H with energy E^^ and 
f{x) = Xj (eS^~^^l^ + 1) is the Fermi function. For ho- 
mogeneous systems (where ^r = 0), the density is 
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where ek = -2Jcos/ca, J\f{E) = {nay^^W^ - E'^)~^' - is 
the density of states, and a is the lattice constant. Note 
that both Eq.([2| and ([3| also apply to free fermions. 

The T = phase diagram of this system is shown 
in Figure 1. There are two QC points, jic = — 2J and 
+2J in the T — ji plane, corresponding to the transi- 
tion from vacuum (n = 0) to a Tomonaga-Luttinger Liq- 
uid (TLL) and from the TLL to Mott phase (n = 1); 
denoted as / and //, respectively. The region T > 
\2J + /i| {QCI in Fig.l) is a quantum critical region, 
where thermal excitations dominate over the excitation 
energy of the system (measured from /i). Similar re- 
gion appears near QC point // by particle- hole symme- 
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FIG. 2: Scaled density (top) and scaled compressibility 
(bottom) vs. /i(x) — II — ^Muj'^x^ at different tempera- 
tures T/J = 0.1,0.15,0.2,0.25 (red, blue, purple, and brown 
curves), for the ID hardcore Bose gas discussed in the text. 
Left and right panels are results for transitions / and //, re- 
spectively. When z is chosen to be 2 and u = 1/2, all scaled 
density (compressibility) intersect at /i//J = —1.98 ( — 1.94) 
and fiiil J = 1.998(1.991). The exact values /ic of these crit- 
ical points are —1 and 1. 
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FIG. 3: Scaled density (top) and scaled compressibility (bot- 
tom) vs. (/i(x) — iic)/T^^ . All the data collapse onto a single 
curve. The parameters are the same as Fig. 2. 



try. Near QCI^ fic = — 2J, and the density and com- 
pressibility scale as \n{ii^T) = Lzi/2(— e^^^^^^^"^) and 
AT^ = Li_i/2(-e(^+2^)/^), where A = h/y/2iTm*T, 
rri' = ?i^/{Jo?) is the effective mass, and Lin{x) = 
X]fe=o 12 ^^/^"^ is the polylog function. 

(C) Demonstration of algorithms in (A): For explicit 
demonstration, we consider a hard core Bose gas with 
N = 265 particles in a trap of frequency u adjusted to 
produce a Mott phase at the trap center, correspond- 
ing to Muj'^a^/{2J) = 0.0001. Choosing z = 2, and 
V = 1/2 (which are the correct values for this model). 
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FIG. 4: Determining /Xc of bulk systems using a finite-size 
scaling scheme: the scaled density and scaled compressibility 
are plotted against /i for different T/J = 0.02,0.03,0.05,0.07 
(red, blue purple, brown, and pink curves) for ID hard core 
bosons in harmonic traps. The value V defined in the text is 
set equal to 0.25. 



the scaled quantities are A^^{x) = rf^{x)T and 
C^^[x) = /^^^(x)Ti/2. In Figure 2, we have plotted 
A^^{x) and C^^{x) versus fi{x) for different tempera- 
tures. We find that for T ^ 0.3J, all the curves at 
different temperatures intersect at /iJ/J = —1.98 and 
J = 1.998, corresponding to the critical points / and 
// in Fig.l. Comparing with the exact value /j^c/J = —2 
and 2, one sees that LDA is accurate to 99% even for a 
system with 10^ particles. At T > 0.3 J, the scaled com- 
pressibility at different temperatures fail to intersect at 
the same point, showing that QC properties appear only 
when T < 0.3J. Similar intersections are found when we 
plot the scaled compressibility with local chemical poten- 
tial (see Figure 2). To construct the scaling function, we 
plot the scaled density and scaled compressibility at dif- 
ferent T against the scaled variable /i(x)/T^/^^^^ around 
QC/ and //. We see from Figure 3 that all data collapse 
onto a single curve near /i*, which is the scaling function 
for these QC points. 

(D) Finite size scaling: Although the small error of 
the LDA scheme (about 1%) is hardly of concern to cur- 
rent experiments, it is a matter of principle whether it 
can be eliminated. This question is not academic, for it 
will also be of experimental relevance as the resolution of 
density measurements continue to improve. Moreover, it 
is also relevant for understanding the results of numer- 
ical calculations. Here, we show that this error can be 
corrected by a (more involved) algorithm that make use 
of finite size scaling rather than LDA. Using renormal- 
ization group arguments J2j [13], it has been shown that 
for a bulk system close to a QC point, if one switches on 
a harmonic potential, then the singular part of a ther- 
modynamic quantity near x = follows certain scaling 



FIG. 5: Scaled density (top) and scaled compressibility (bot- 
tom) vs. {fi — iic)/T^^^^ with vz — 1. All curves in Figure 4 
collapse onto a single curve. The parameters are the same as 
Fig. 4. 

relations. The validity of these scaling relations has been 
well tested for quantum spin systems [121 US]. For ex- 
ample, under a scale change 6, the singular part of the 
density below dc scales as 

n,(/x,T,w2^x) = b-^''+'^+^/''g{fibi,T¥,Lo%y,x/b), 

(4) 

where ji = fi — (ic^ ^ = 2 + 1/i^, and g is scaling function. 
Choosing Th^ = 1, the scaled density A{jj.^T] D^x.) = 
T-i-^/^+i/^^n5(T,/i,cc;^x) at point x satisfies 

A(/i, T, x) = ^ (/i/T^ , uj^t'^ ,xTiy (5) 

Defining D = uo'^ /T^ , Eq.([5| becomes 

A(/i,T|Ax) =^ (/i/Ti, AxT^) , (6) 

where ]4(/i,T|D,x) = A{/j.,T, DT^^ We then have 
A{/j.,T\D,0) = g (^fl/T^,D,Oy Hence, if we plot A ver- 
sus /i at X = for different T with D held fixed, then dif- 
ferent curves with different T (but same D) will intersect 
at jiic since A{jiic^ T\D^ 0) = ^ (0? D^O). A similar analysis 
shows scaled compressibility T~^~^/^+^/^^Ai:5(T, /i, cj^, x) 
at constant D also intersect at /ic- 

For ID hard core bosons, z = 2, = 1/2, = 
n, we have A(/i,T|D,x) = n'''' {ii,uj,T]x)T-^^'^ near 
X = 0. Our scheme is to plot n^^ {/j.,uj,T;0)T~^^'^ ver- 
sus /i for a family of samples with different T and uj 
but with identical values for the ratio D = uj'^/T^. 
These plots, together with similar plots for the com- 
pressibility, are given in Figure 4, with the ratio V = 
{\Muj'^a^/J)/{T/J)yl^ = 0.25. We find that, to machine 
accuracy, different curves intersect exactly at jHc/ J = 2.0. 
Moreover, if we plot rf^{ii^ T; cj, 0)T~^/^ against the vari- 
able u = {ji — iic)/T^^^^ ^ we see in Fig. 5 that all curves 
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collapse to the same curve, which is the scaling function 
g{u,D,0). 

By holding D fixed, we need to know the densities 
{n^^(/i, T; uj)} for a range of T, /i and uj values. Thus, this 
scheme is clearly more laborious than the LDA scheme. 
Still, it is of conceptual importance to demonstrate that 
the /ic of bulk systems can be obtained exactly without 
errors. The recent work by Salomon's group [9 mapping 
out the equation of state of a unitary gas indicates that 
the construction of the entire family {n^^(//, T; uj)} is fea- 
sible. Our scheme will also be useful in numerical studies 
of trapped atoms, especially for unsolved problems where 
the existence of quantum criticality is in question. 

(E) Finite temperature phase transition: Previously, 
we [7] have pointed out that within LDA, the phase 
boundary fic{T) of a finite temperature continuous phase 
transition will show up as a kink in the in the compress- 
ibility n^^{x) = dn^^ {x) / d ii{x) at position x* such that 
IJic(T) = ii{x*). As with the quantum critical case dis- 
cussed in main text, LDA will contain a systematic error 
due to finite size effects, even though it is very small. 
This error can also be corrected in a similar way. The 
analog of Eq. Q for a finite temperature phase transition 
is[l6] 

n,(T, /i, a;^ x) = b-^^'^^ns{-p{T)b'^\oJ%^^/b), (7) 

where /i(T) = ii—iic{T) and v is the correlation length ex- 
ponent for finite T transition. (For example, v = 0.67 for 
the 3D xy universality class.) Setting = 1, Eq. Q 
implies the scaling form 

n,(T,/i,cJ^x) =cj^^(cj-2/^"/i(T),a;2/^x), (8) 

where x = {2/y){d — l/u) and is a scaling 
function. Similar analysis gives ^^^(T, /i, cj^, x) = 
uj^'^^y^^^-'^^^^T{uj-'^^y''Jl{T),u'^^yx.). The fuh density n, 
however, contains both singular (n^) and regular (n^) 
contributions; i.e. n = + ng. The unknown regu- 
lar contribution can be eliminated by measuring the 
difference in density profiles at different trap frequencies 
UJ uji relative to some reference value uji. Fixing cji, 
the function 

^ n(^,M,.^0)-n(^,M,.^0) 

will assume the value J^(0, 0) at /i = Mc(^) for all uj. This 
means that if we plot Bt{/J^;cjj) as a function of /i, then 
different curves for different cu will intersect at the exact 
phase boundary /ic(^)- 

The fact that the entire phenomenology of quantum 
criticality of bulk systems - its existence, its dynamical 
critical exponent, its scaling function, and the location of 
the quantum critical point - can all be deduced from the 
T 7^ density profile of trapped gases as demonstrated 
here further adds to the list of valuable information de- 
ducible from density measurements At present, the 



field of quantum gases is moving rapidly in the direction 
of high precision measurements, a direction essential for 
realizing the full power of quantum simulation [7]. The 
extraordinarily high resolution for density imaging re- 
cently achieved by Greiner's group [HI [15] shows that 
the accuracy needed for deducing bulk properties from 
trapped gases has been achieved. The recent success 
of Salomon's group in deducing the phase diagram of 
a unitary gas [8, 9 demonstrates the feasibility of our 
algorithms. All these developments strongly suggest the 
power of quantum simulation is ready to be harvested. 
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